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^  A  FORTRAN  COntrol  Program  for  Engineering  Synthesis 
(COPES)  has  been  developed  for  solving  engineering  design 
problems.  The  program  maximizes  or  minimizes  a  numerically 
defined  objective  function  subject  to  a  set  of  inequality 
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tions  for  constrained  function  minimization.  Additionally, 
approximation  techniques  are  available  for  use  in  optimiza¬ 
tion,  and  trade-off  studies  may  be  performed.  A  simple 
design  example  demonstrates  the  program  capabilities. 
Programming  guidelines  are  presented  followed  by  sample 
input  data  and  output  for  each  program  option. 
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I.  INTRODUCTION 


Most  design  processes  require  the  minimization  or  maximi¬ 
zation  of  some  parameter  which  may  be  called  the  design  ob¬ 
jective.  For  the  design  to  be  acceptable,  it  must  satisfy  a 
variety  of  physical,  aesthetic,  economic  and,  on  occasion, 
political  limitations  which  are  referred  to  here  as  design 
constraints.  While  part  of  the  design  problem  may  not  be 
easily  quantified,  most  of  the  design  criteria  can  be  des¬ 
cribed  in  numerical  terms. 

To  the  extend  that  the  problem  can  be  stated  in  numerical 
terms,  a  computer  program  can  be  written  to  perform  the 
necessary  calculations.  For  this  reason,  computer  analysis 
is  commonplace  in  most  engineering  organizations.  For  example 
in  structural  design,  the  configuration,  materials  and  loads 
may  be  defined  and  a  finite  element  analysis  computer  code 
is  used  to  calculate  stresses,  deflections  and  other  response 
quantities  of  interest.  If  any  of  these  parameters  are  not 
within  the  prescribed  bounds,  the  engineer  may  change  the 
structural  member  sizes  and  re-run  the  program.  The  computer 
code  therefore  provides  only  the  analysis  of  a  proposed  design 
with  the  engineer  making  the  actual  design  decisions.  This 
approach  to  design,  which  may  be  called  computer-aided  design, 
is  commonly  used  today. 


Another  common  use  of  analysis  codes  is  in  trade-off 
studies.  For  example,  an  aircraft  trajectory  analysis  code 
may  be  run  repetitively  for  several  payloads,  calculating  the 
aircraft  range,  to  determine  the  range-payload  sensitivity. 

A  logical  extension  to  computer-aided  design  is  fully 
automated  design,  where  the  computer  makes  the  actual  design 
decisions,  or  to  perform  trade-off  studies  with  a  minimum  of 
man-machine  interaction.  The  purpose  of  the  COPES  program 
is  to  provide  this  automated  design  and  trade-off  capability. 
The  user  must  provide  a  FORTRAN  analysis  program  for  analysis 
of  the  particular  problem  being  considered.  This  analysis 
program  is  written  according  to  a  simple  set  of  guidelines 
so  that  it  can  be  easily  coupled  to  the  COPES  program  for 
automated  design  synthesis. 

This  document  describes  the  capabilities  of  the  COPES 
program  and  its  usage.  A  simple  design  example  is  first  pre¬ 
sented  to  demonstrate  the  program  capabilities.  Guidelines 
are  given  for  writing  analysis  codes  which  can  be  coupled 
directly  to  COPES.  Finally  the  data  organization  is  outlined 
and  sample  data  is  presented. 

The  most  commonly  used  option  of  the  COPES  program  will 
be  for  design  optimization.  Two  approaches  are  available 
for  this  purpose.  The  first  (NCALC  =  2)  is  direct  optimiza¬ 
tion  of  the  function  by  the  CONMIN  optimization  sub-program 
[Ref.  1] .  An  alternative  to  this  is  through  the  use  of 
approximation  techniques  (NCALC  *  6,  Ref.  2).  This  second 


option  is  usually  more  efficient  for  problems  of  under  6 
design  variables  but  which  requires  costly  analysis,  espe¬ 
cially  when  multiple  optimizations  are  to  be  performed. 


DESIGN  EXAMPLE 


Assume  it  is  required  to 
design  the  cantilevered  beam 
shown  in  Fig.  1.  The  objec¬ 
tive  is  to  find  the  minimum 
volume  of  material  which  will 
support  the  concentrated  load. 
That  is; 


P*1 0,000  lb. 

! 


O' 


L*200  in. 


I 


E*30xl0°  psi 


Minimize  Volume  =  B«H«L 


CD 


Figure  1  -  Cantilevered 
Beam 


The  bending  stress  in  the  beam  must  not  exceed  20,000  psi; 
-  _  Me  _  6PL. 


TH2  1  20  ,000 


(2) 


The  shear  stress  must  not  exceed  10,000  psi; 


v  "  1  I  "  IlH  1  10’000 


(3) 


and  the  deflection  under  the  load  must  not  exceed  one  inch; 

(4) 


.  .  PL’  .  4 PL 3  .  „ 

3lT  EOT*  1  1-0 


Additionally,  geometric  limits  are  imposed  on  the  beam  size 
so  that; 


0.S  <  B  <  5.0  (5) 

1.0  <  H  <  20.0  (6) 

H/B  <  10.0  (7) 


Equation  (1),  the  design  objective,  provides  a  measure  of 
the  efficiency  of  the  design,  while  Eqs .  (2-7)  define  the 
criteria  which  the  structure  must  satisfy,  and  are  referred 
to  as  constraints. 

This  design  problem  may  be  stated  in  general  form; 

Minimize  F(X) 

Subject  to: 

Gj  (X)  £  0  i  *  1  ,m 

xi  i  xi  1  xi 

where  X  is  a  vector  containing  the  design  variables  B  and  H 
and  G . (X)  are  the  constraints  defined  by  Equations  (2-7). 
There  are  eight  constraints  on  the  design.  The  objective 
and  constraints  are  functions  of  the  design  variables. 

While  this  problem  is  straightforward,  its  solution  is 
not  trivial  because  it  is  not  known  which  constraints  will 
be  critical  (i.e.  which  G ■ (X)  -  0)  at  the  optimum. 

One  design  approach  is  to  assume  that  the  bending  and 
the  displacement  constraints  are  critical,  so  that  Equations 
(2)  and  (4)  are  equalities,  and  then  solve  the  two  simulta¬ 
neous  equations  for  the  design  variables  B  and  H.  The 
remaining  constraints  are  then  checked  to  be  sure  they  are 
satisfied,  and  if  not,  they  are  used  to  obtain  a  new  design 
satisfying  all  constraints.  When  a  design  is  obtained  where 
two  Gj (X)  *  0  and  the  remaining  Gj  (X)  £  0,  it  is  an  accept¬ 
able  design. 
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The  design  obtained  here  yields  a  volume  of  6750  and 
satisfies  all  constraints.  Note,  however,  that  the  objec¬ 
tive  function  of  Equation  (1)  played  no  part  in  the  choice 
of  design  variables  so  that,  using  this  approach,  there  is 
no  assurance  that  a  minimum  volume  design  will  be  produced 
[Ref.  3].  More  importantly,  it  is  desirable  to  devise 
techniques  for  optimum  design  of  systems  which  may  be  de¬ 
fined  by  more  than  two  variables  and  by  much  more  complex 
analysis,  where  the  optimum  design  cannot  be  determined  by 
inspection. 

The  COPES  program  provides  this  general  capability  by 
the  use  of  the  optimization  program  CONMIN  [Ref.  1].  To  use 
this  design  capability,  a  FORTRAN  code  must  be  provided 
which  will  calculate  the  various  parameters.  In  writing  the 
analysis  code;  1)  it  is  written  in  subroutine  form  with 
SUBROUTINE  ANALIZ  (ICALC)  as  the  main  routine,  2)  it  is  seg¬ 
mented  into  INPUT,  EXECUTION  AND  OUTPUT  and  3)  all  parameters 
which  may  be  design  variables,  object  functions  or  con¬ 
straints  are  contained  in  a  single  labeled  common  block  called 
GLOB CM. 

To  demonstrate  the  simplicity  with  which  a  design- 
oriented  analysis  code  can  be  written,  the  following  FORTRAN 
subroutine  was  produced  for  the  analysis  of  the  cantilevered 
beam  in  Figure  1. 


SUBROUTINE  ANAL I Z  (I CALC) 

COMMON  /GLOB CM/  B ,H , VOL , BSTRES ,SHRSTR, DELTA ,HB ,E ,AL 
IF  (ICALC.GT.l)  GO  TO  10 
C  ---  INPUT  OR  INITIALIZATION. 

B-2 . 5 
H=10 . 

P-10,000. 

E-30.E+6 
AL-200 . 

WRITE  (6 , 30) AL , P ,E  ,B  ,H 
RETURN 

C  ---  EXECUTION. 

10  CONTINUE 

VOL=AL*B*H 

BSTRES-6 . *P*AL/ (B*H*H) 

SHRSTR-1 . 5*P/ (B*H) 

DELTA-4 . *P* (AL**3) /(E*B*(H**3)) 

HB-H/B 

IF  (ICALC.LT. 3)  RETURN 
C  ---  PRINT  RESULTS. 

20  WRITE  (6,30) AL,P,E,B,H 

WRITE  (6,40) VOL , BSTRES ,SHRSTR, DELTA, HB 
30  FORMAT (// ///5X, 17HCANTILEVERED  BEAM//5X, 8HAL  =,F9.3/5X, 

*  8HP  = ,E12 . 5/5X, 8HE  =,E12 .5//5X,8HB  =,F9.3/5X, 

*  8HH  = ,F9 . 3) 

40  FORMAT (/5X, 8HVOL  -,F9 . 3//5X, 8HBSTRES- ,E12 . 5/5X, 
8HSHRSTR- 

*  E12 . 5/5X, 8HDELTA- ,E12 . 5/5X, 8HH/B  =,F9.3) 

RETURN 

END 

This  routine  may  be  executed  as  a  simple  analysis  program 
by  the  following  main  program; 

C  MAIN  PROGRAM  TO  EXECUTE  SUBROUTINE  ANALIZ. 

DO  10  I  -  1,3 
10  CALL  ANALIZ  (I) 

STOP 

END 

Moreover,  ANALIZ  can  be  coupled  directly  to  the  COPES 
program  to  perform  this  same  function,  or  to  perform  optimiza¬ 
tion  or  trade-off  studies. 

This  subroutine  was  coupled  to  COPES  and  an  optimization 
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was  performed  (NCALC  =  2)  to  yield  the  following  design; 

CANTILEVERED  BEAM 

L  *  200.00 

P  =  .  10000E+05 

E  *  . 30000E+0  8 

B  =  1.82 

H  =  18.16 

VOL  *  6607.61 

BSTRES  =  .  20002E+05 

SHRSTR  *  . 45402E  +  03 

DELTA  =  . 97906E+00 

H/B  =9.98 

The  parameters  L,  P,  E,  B  and  H  were  input  to  the  pro¬ 
gram  and  the  remaining  parameters  were  calculated.  The  design 
variables  B  and  H  were  changed  during  the  optimization  process 
to  obtain  this  design. 

This  design  was  achieved  with  45  calls  to  ANALIZ  with 
ICALC  =  2  (45  analyses).  This  design  could  surely  have  been 
found  with  fewer  analyses  were  it  performed  by  hand  calcula¬ 
tions.  However,  once  having  written  the  analysis  subroutine, 
numerous  other  designs  may  be  obtained  for  different  materials, 
loading,  or  design  stresses  with  only  minimal  effort. 
Furthermore,  the  design  obtained  here  of  volume  *  6607.6  is 
very  near  the  theoretical  optimum  of  6603.9,  while  that  which 
satisfied  stress  and  displacement  constraints  simultaneously 
was  not. 

This  very  simple  design  example  underscores  the  power  of 
numerical  optimization  techniques  and  the  ease  with  which  they 
may  be  applied.  The  key  to  efficient  use  of  the  COPES  program 
is  the  requirement  that  the  ANALIZ  code  be  written  in  a  standard 
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format  so  that  it  may  be  coupled  to  COPES  without  modifica¬ 
tion.  The  following  section  contains  guidelines  for  writing 
a  design-oriented  analysis  code. 


II.  PROGRAMMING  GUIDELINES 


In  developing  any  computer  code  for  engineering  analysis, 
it  is  prudent  to  write  the  code  in  such  a  way  that  it  is 
easily  coupled  to  a  general  synthesis  program  such  as  COPES. 
Therefore,  a  general  programming  practice  is  outlined  here 
which  in  no  way  inhibits  the  use  of  the  computer  program  in 
its  traditional  role  as  an  analytic  tool,  but  allows  for  sim¬ 
ple  adaptation  to  COPES.  This  approach  is  considered  good 
programming  practice  and  provides  considerable  flexibility 
of  design  options.  Only  five  basic  rules  must  be  followed: 

I.  Write  the  code  in  subroutine  form  with  the  primary 
routine  called  as;  SUBROUTINE  ANALIZ  (ICALC)  The 
name  ANALIZ  is  compatable  with  the  COPES  program  and 
ICALC  is  a  calculation  control.  Note  that  subroutine 
ANALIZ  may  call  numerous  other  subroutines  as  re¬ 
quired  to  perform  the  necessary  calculations. 

II.  Segment  the  program  into  INPUT,  EXECUTION  and  OUTPUT. 
The  calculation  control,  ICALC,  will  determine  the 
portion  of  the  analysis  code  to  be  executed. 

ICALC  ■  1;  the  program  reads  all  data  required  to  perform 
the  analysis.  Also,  any  initialization  of  constants  which 
will  be  used  repetitively  during  execution  is  done  here. 

This  initial  input  information  is  printed  here  for  later  re¬ 
ference  and  for  program  debugging. 


I CALC  *2:  the  program  performs  the  execution  phase  of 
the  analysis  task.  No  data  reading  or  printing  is  done  here, 
except  on  user-defined  scratch  disc.  Data  may  be  printed 
here  during  program  debugging,  in  which  case  it  should  be 
controlled  by  a  print  control  parameter  which  is  read  during 
input.  In  this  way,  this  print  may  be  turned  off  after  the 
program  is  debugged,  but  may  be  used  again  during  future 
program  expansion  and  debugging.  The  reason  that  printing 
is  not  normally  allowed  during  execution  is  that  when 
optimization  is  being  done,  the  code  will  be  called  many 
times  with  ICALC  -  2,  resulting  in  voluminous  print. 

ICALC  *  3:  the  results  of  the  analysis  are  printed. 

Also  the  essential  input  parameters  which  may  have  been 
changed  during  optimization  should  be  printed  here  for  easy 
reference.  In  some  cases,  so  much  information  is  generated 
when  ICALC  *  2  that  it  would  be  inefficient  to  store  it 
internally  or  on  disc  for  printing  when  ICALC  =3.  In  this 
case,  it  may  be  desirable  to  actually  execute  when  ICALC  *  3 
with  a  print  code  turned  on  to  print  his  intermediate  infor¬ 
mation.  However,  this  approach  should  be  avoided  because 
it  requires  an  additional  complete  execution  of  the  program. 

In  summary,  when; 

ICALC  *  1  Read  input  data. 

ICALC  *  2  Execute  the  analysis. 

ICALC  ■  3  Print  the  results. 


III.  Store  all  parameters  which  may  be  design  variables, 
objective  functions  or  constraints  in  a  single 
labeled  common  block  called  GLOBCM.  The  order  in 
which  they  are  stored  is  arbitrary.  A  listing  of 
the  COPES  program  should  be  checked  to  see  how  many 
parameters  my  be  stored  in  GLOBCM  (the  dimension  of 
ARRAY) .  Initial  distribution  of  COPES  allows  for 
1500  parameters. 

IV.  During  execution  or  output,  no  parameters  which  are 
read  during  input  should  be  updated.  For  example, 
if  variable  X  is  initialized  during  input,  the 
execution  segment  must  not  update  X  such  as  X  =  X  + 
3.2.  Instead  a  new  variable,  Y  *  X  +  3.2  should  be 
defined. 

V.  Write  all  programs  in  standard  language,  avoiding 
machine  dependent  capabilities  such  as  multiple 
entry  point  (IBM) ,  DEFINE  statements  (UNIVAC)  and 
seven  letter  FORTRAN  names  (CDC) .  While  this  guide¬ 
line  is  not  essential  to  the  use  of  the  analysis 
code  within  the  COPES  program,  it  makes  the  analysis 
code  much  more  transportable  between  different 
computer  systems,  a  capability  which  easily  justi¬ 
fies  a  slight  reduction  in  efficiency  on  a  given 
machine. 

Adherence  to  these  guidelines  not  only  leads  to  a  more 
readable  and  machine  independent  computer  code,  but  allows 
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this  code  to  be  coupled  to  the  COPES  program  without 
modification. 

Having  written  the  analysis  code,  it  may  be  executed 
either  with  a  simple  main  program  or  within  the  COPES  program 
to  perform  the  analysis.  To  insure  that  guideline  IV  is 
followed,  the  following  main  test  program  is  recommended. 

Note  that  this  program  calls  ANNALIZ  twice  with  ICALC  *  2  and 
ICALC  *  3,  to  show  that  the  same  result  is  obtained  repetitively 

C  MAIN  PROGRAM  TO  CHECK  SUBROUTINE  ANALIZ, 

C  READ,  EXECUTE  AND  PRINT. 

DO  10  ICALC  -  1,  3 
10  CALL  ANALIZ  (ICALC) 

C  EXECUTE  AND  PRINT  AGAIN  TO  BE  SURE  THE  RESULTS 

C  DO  NOT  CHANGE. 

DO  20  ICALC  -  2,  3 
20  CALL  ANALIZ  (ICALC) 

STOP 

END 

This  program  was  executed  with  the  ANALIZ  subroutine  for 

the  cantilevered  beam  example  to  yield  the  following  result: 

CANTILEVERED  BEAM 

AL  -  200.00 

P  -  . 10000E+05 

E  -  . 30000E+08 


CANTILEVERED  BEAM 


AL  =  200.00 

P  *  . 10000E+05 

E  =  . 30000E+08 

B  *  2.00 

H  =  5.00 

VOL  *  2000.00 

BSTRES  «  . 24000E+06 

SHRSTR  *  . 15000E+04 

DELTA  =  . 42667E+02 

H/B  *  2.50 

CANTILEVERED  BEAM 


AL 

P 

E 

B 

H 

VOL 

BSTRES 

SHRSTR 

DELTA 

H/B 


*  200.00 
=  . 10000E+05 
=  . 30000E+08 
=  2.00 
»  5.00 
“  2000.00 
*  . 24000E+06 
=  . 15000E+04 
»  . 42667E+02 
=  2.50 


This  design  was  used  as  the  initial  design  for  the  optimi¬ 
zation  presented  previously.  Note  that,  while  the  volume  here 
is  lower  than  the  optimum,  the  bending  stress  and  displacement 
each  exceed  the  imposed  limits  by  more  than  an  order  of 
magnitude. 

Once  the  analysis  code  has  been  written,  it  can  be  coupled 
to  the  COPES  program  without  modification.  If  it  is  desired 
to  perform  a  simple  analysis  using  COPES,  only  three  data 
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cards  are  required  for  the  COPES  program,  namely  a  TITLE  card, 
a  control  parameter,  NCALC  *  1,  and  an  END  card.  If  the 
optimization  or  parametric  analysis  (sensitivity)  capabilities 
of  COPES  are  to  be  used,  additional  data  must  be  read.  This 
data  will  identify  which  parameters  in  the  global  common  block, 
GLOBCM,  are  used.  To  set  up  the  COPES  data,  the  user  must 
have  a  basic  understanding  of  how  the  data  in  the  global  com¬ 
mon  block  is  accessed  by  COPES.  This  is  outlined  in  the 
following  section. 


III.  DATA  MANAGEMENT 


In  order  to  perform  design  operations,  the  COPES  program 
must  access  the  data  in  common  block  GLOBCM.  This  is  done  by 
defining  the  location  in  GLOBCM  where  a  specified  parameter 
resides.  For  example,  consider  the  common  block  for  the  can¬ 
tilevered  beam; 

COMMON/ GLOBCM/B,H, VOL, BSTRES .SHRSTR, DELTA, HB,E,AL 

The  volume  of  material,  VOL,  is  the  third  parameter  in 
the  common  block;  that  is,  it  resides  in  location  3,  referred 
to  as  the  global  location  number.  Similarly  the  bending 
stress,  BSTRES,  is  in  global  location  4  and  the  beam  width  is 
in  global  location  1.  Thus,  the  parameters  are  referred  to 
by  their  respective  location  numbers  in  global  common. 

For  convenience  in  preparing  data  for  the  COPES  program, 
a  simple  "CATALOG"  of  parameters  may  be  defined.  For  the 
cantilevered  beam,  this  catalog  would  be; 


GLOBAL 

LOCATION 

FORTRAN 

NAME 

DEFINITION 

1 

B 

Beam  width 

2 

H 

Beam  height 

3 

VOL 

Volume  of  Material 

4 

BSTRES 

Maximum  bending  stress 

5 

SHRSTR 

Maximum  shear  stress 

6 

DELTA 

Deflection  under  the  load 

7 

HB 

Ratio,  H/B 

20 


GLOBAL 

LOCATION 


FORTRAN 

NAME 


DEFINITION 


I 


8 

E 

Young’s  modulus 

9 

AL 

Length  of  beam 

As 

another  example,  consider  a  global  common  block  con 

taining 

arrays ; 

COMMON/ GLOBCM/A,  Y(10),  Q 

,  C(2 ,2)  ,  H 

The 

variable  catalog  for 

this  common  block  is; 

GLOBAL 

FORTRAN 

LOCATION 

NAME 

DEFINITION 

1 

A 

Area 

2 

Y(10) 

Vector  of  y-coordinates 

12 

Q 

• 

13 

C(2 ,2) 

• 

17 

H 

etc. 

The  dimensions  are  given  with  the  FORTRAN  name  as  a 
reminder  that  the  parameter  is  an  array.  In  this  case,  the 
third  parameter  in  the  Y  array  is  in  global  location  4. 
Remembering  that  arrays  are  stored  column  by  column  the 
C(l,2)  array  location  is  in  global  location  15. 

It  will  be  seen  that  identifying  parameters  according 
to  their  location  in  GLOBCM  provides  a  great  deal  of  flexi¬ 
bility  in  using  the  COPES  program  for  design. 

Appendix  A  contains  blank  forms  for  writing  the  variabl 


catalog  for  the  user's  particular  problem. 


IV.  COPES  TERMINOLOGY 


The  copes  program  currently  provides  six  specific 
capabilities ; 

1.  Simple  analysis,  just  as  if  COPES  was  not  used. 

2.  Optimisation  -  Minimization  or  maximization  of  one 
calculated  function  with  limits  imposed  on  other 
functions  using  the  CONMIN  subprogram. 

3.  Sensitivity  analysis  -  the  effect  of  changing  one 
or  more  design  variables  on  one  or  more  calculated 
functions . 

4.  Two- variable  function  space  -  analysis  for  all 
specified  combinations  of  two  design  variables. 

5.  Optimum  sensitivity  -  same  as  sensitivity  analysis 
except  at  each  step,  the  design  optimized  with 
respect  to  the  remaining  independent  design  vari¬ 
ables  . 

6.  Approximate  optimization  -  optimization  using 
approximation  techniques.  Usually  more  efficient 
than  standard  optimization  for  up  to  6  design 
variables  or  if  multiple  optimizations  are  to  be 
performed. 

In  defining  the  data  required  to  execute  the  COPES  pro¬ 
gram,  the  following  definitions  are  useful: 

Design  Variables  -  Those  parameters  which  the  optimiza¬ 
tion  program  is  allowed  to  change  in  order  to  improve  the 
design.  Design  variables  appear  only  on  the  right  hand  side 


of  equations  in  the  analysis  program.  COPES  considers  two 
types  of  design  variables,  independent  and  dependent.  If 
two  or  more  variables  are  always  required  to  have  the  same 
value  or  be  in  a  constant  ratio,  one  is  the  independent 
variable  while  the  remaining  are  dependent  variables.  For 
example,  if  the  height  is  required  to  be  10  times  the  width 
of  the  cantilever  beam,  B  would  be  the  independent  variable 
while  H  would  be  the  dependent  variable. 

Objective  Function  -  The  parameter  which  is  to  be  mini¬ 
mized  or  maximized  during  optimization.  Also  the  parameters 
calculated  as  functions  of  specified  design  variables  during 
a  sensitivity  or  two-variable  function  space  study.  Objective 
functions  always  occur  on  the  left  side  of  equations,  unless 
the  objective  function  is  also  a  design  variable  (the  beam 
height  may  be  minimized  as  an  objective  function  if  it  is 
also  a  design  variable.  In  this  case  the  minimum  height  is 
found  for  which  no  constraints  are  violated).  An  objective 
function  may  be  linear  or  non-linear,  implicit  or  explicit, 
but  must  be  a  function  of  the  design  variables  to  be  meaning¬ 
ful. 

Constraint  -  Any  parameter  which  must  not  exceed  speci¬ 
fied  bounds  for  the  design  to  be  acceptable.  Constraint 
functions,  always  appear  on  the  left  side  of  equations.  Just 
as  for  objective  functions,  constraints  may  be  linear  or 
non-linear,  implicit  or  explicit,  but  must  be  functions  of 
the  design  variables. 
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Constraint  Set  -  A  group  of  constraints  which  appear 


consecutively  in  the  global  common  block  and  which  all  have 
the  same  limits  imposed.  This  is  a  convenience  which  allows 
several  constraints  to  be  identified  with  a  minimum  of  data. 

Global  Common  -  Common  block  GLOBCM  containing  design 
information. 

Global  Location  -  Location  of  a  particular  parameter  in 
GLOBCM. 
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V.  COPES  DATA 


The  COPES  program  reads  data  from  unit  5  and  writes  output 
on  unit  6.  Units  20  and  40  are  used  as  scratch  files.  The 
scratch  file  numbers  may  be  changed  by  changing  two  cards  at 
the  beginning  of  the  COPES  program. 

In  order  to  execute  the  COPES  program  it  is  necessary  to 
provide  formatted  data  for  COPES,  followed  by  data  for  the 
ANALIZ  program  which  is  coupled  to  COPES.  This  section  de¬ 
fines  the  data  which  is  required  by  COPES.  The  data  is  seg¬ 
mented  into  "BLOCKS"  for  convenience.  All  formats  are  alpha¬ 
numeric  for  TITLE,  and  END  cards,  F10  for  real  data  and  110 
for  integer  data. 

The  COPES  data  begins  with  a  TITLE  card  and  ends  with  an 
END  card.  This  is  followed  by  data  to  be  read  by  the  user 
supplied  subroutine  ANALIZ  or  when  ICALC  =  1. 

Comment  cards  may  be  inserted  anywhere  in  the  COPES  data 
stack  prior  to  the  END  card,  and  are  identified  by  a  dollar 
sign  ($)  in  column  1. 

Data  coding  forms  are  provided  in  Appendix  B. 


UNFORMATTED  DATA  INPUT 


While  the  user's  sheet  defines  COPES  data  in  formatted 
fields  of  ten,  the  data  may  actually  be  read  in  a  simplified 
fashion  by  separating  data  by  commas  or  one  or  more  blanks. 

If  more  than  one  number  is  contained  on  an  unformatted  data 
card,  a  comma  must  appear  somewhere  on  the  card.  If  expo¬ 
nential  numbers  such  as  2.5+10  are  read  on  an  unformatted  card 
there  must  be  no  embedded  blanks.  Unformatted  cards  may  be 
intermingled  with  formatted  cards.  Real  numbers  on  an  un¬ 
formatted  card  should  have  a  decimal  point. 

Examples : 


Unformatted  data; 

5,7,1. 3,1. 0+20,0, -5.1 
5,7,1.3,1.0+20, ,-5.1 
5  7  1.3  1.0+20,,  -5.1 
5  7  1.3,  1.0+20  0  -5.1 

Equivalent  formatted  data; 
col-+  10  20  30 


Equivalent 

formatted 

data; 

col-*- 

10 

20 

30 

2 

Note: 

This  data  has  been 

right 

2 

3 

justified. 

2  3 

Note : 

This  data  contains 

no  commas , 

so  it  is  assumed  to 
ted  already. 

be  format 

Unformatted 

data; 

1,2,3 

,4,5,6 

,7,8,9,10, 

11 

Equivalent 

formatted 

data; 

col-*- 

10 

20 

30 

40  50  60 

70 

Note  that  two  formatted  data  cards  are  created  here 
Unformatted  data; 

1,2, 3, 4,5, 6, 

7,8,9,10,11 

Equivalent  formatted  data; 

col-*  10  20  30  40  50  60 


Note  that  the  above  two  examples  do  not  produce  the  same  for 
matted  data  cards. 


DESCRIPTION:  Program  Control  Parameters 
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Two  variable  function  space.  Data  of  blocks  A,  B,  and  R-V  is  required 
Remaining  data  is  optional. 


cont.  NCALC 


DATA  BLOCK  C  OMIT  IF  NDV  -  0  IN  BLOCK  B 


Print  above  plus  each  proposed  design  vector,  objective  function 
and  constraint  values  during  the  one -dimensional  search. 


DESCRIPTION:  Floating  point  optimization  program  parameters. 
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DATA  BLOCK  G  OMIT  IF  NDV  =  0  IN  BLOCK  B 


DATA  BLOCK  H  OMIT  IF  NDV  =  0  IN  BLOCK  B 


LCON:  Linear  constraint  identifier  for  this  constraint  set 

LCON  =  1  indicates  linear  constraints. 


BL:  Lower  bound  on  the  constrained  variables.  If  BL.LT. -1 .OE+15 ,  no 
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If  NXS  =  NXFS  =  0,  NXA  is  defaulted  to  NXA  =  1,  even  if  it  is  read  as  zero. 
Also,  a  second  vector  of  design  variables  is  automatically  defined  by  COPES 
to  yield  two  independent  designs  to  start  the  optimization. 
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DATA  BLOCK  L  OMIT  IF  NXAPRX  8  0  IN  BLOCK  B  OR  INXLOC  -  0  IN  BLOCK  J 


FNXAPRX  »  0  IN  BLOCK  B  OR  INFLOC 


E)  followed  by  the  global  locations  of  the  constrained  parameters  (ICON-JCON 
in  BLOCK  I) . 
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DATA  BLOCK  R  OMIT  IF  N2VAR  =  0  IN  BLOCK  B 


DATA  BLOCK  T  OMIT  IF  N2VAR  =  0  IN  BLOCK  B 


DATA  BLOCK  U  OMIT  IF  N2VAR  *  0  IN  BLOCK  B 


VII.  SAMPLE  COMPUTER  SOLUTIONS 


Sample  solutions  are  presented  here  for  the  various  COPES 
program  options  defined  by  the  parameter  NCALC  in  DATA  BLOCK 
B.  The  ANALIZ  routine  given  in  the  design  example  section  was 
used  to  produce  the  results  given  here.  Note  that  SUBROUTINE 
ANALYZ  does  not  read  data,  so  only  COPES  data  is  required  here 
In  the  usual  case  where  data  is  read  as  input  to  the  analysis 
routine,  this  data  would  follow  directly  after  the  COPES  'END' 
card. 

The  output  from  a  COPES  program  execution  includes  a  title 
page  followed  by  a  copy  of  the  input  data.  Then  the  required 
program  executions  are  performed  and  final  output  information 
is  printed.  Figures  2-7  contain  output  for  all  options  of 
COPES  as  follows: 


NCALC  FIGURE 
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PROBLEM  SOLVED 

Analysis  only 

Optimization  by  CONMIN 

Sensitivity  study 

Two-variable  function  space  study 

Optimum  sensitivity 

Optimization  using  approximation 
techniques 


Note  that  data  for  the  COPES  program  is  read  depending  on 
the  value  of  the  parameters  NDV,  NSV,  N2VAR  and  NXAPRX  in  DATA 
BLOCK  B.  The  actual  program  execution  is  determined  by  the 
value  of  NCALC.  Therefore,  data  may  be  read  for  all  program 
options  even  though  only  one  option  of  the  program  is  to  be 
executed.  Figure  8  is  a  copy  of  the  combined  data  for  all  the 
previous  examples.  This  data  may  be  used  for  any  program  op¬ 
tion  simply  by  changing  the  value  of  NCALC.  Figure  9  is  a 
copy  of  this  same  data  using  the  simplified  data  input  mode 
and  Figure  10  is  the  same  data  without  the  comment  cards. 


56 


C A  l/M/l(/)t/) 

CO  (/>  to 

CO  CO  CO 

CO  CO  CO 

(/I  CO  CO 

CO  JO  CO 

CO  COCO  CO  CO 


CO 


CO 


u> 

Ui 

IU 


UI 

Ui 

Ui 

UJ 


z 

o 

CO 


Ui  Ui  Ui 

o 

z 

o 

UiUiUJuiUiUiUi 

o 

► 

o 

z 

fl£ 

CO 

Ui 

4 

a.  a  a  a. 

a.  a 

a 

ae 

•J 

to 

K 

a  “• 

M 

4 

a.  ct 

o 

a 

t- 

CO 

a 

u  <* 

* 

a.  u. 

«J 

u- 

z 

•m* 

mi 

_j 

0.0- 0.0.  a  a.  a 

4 

UI 

VI 

M 

z 

a 

«* 

i- 

ec 

a 

z 

UOUOdOU 

* 

a 

a  u 

t- 

UJ 

4 

Ui 

U  UI 

Ui 

u  u 

z 

Ui 

Ui 

li¬ 

u  u 

es 

Ui  u 

Ui 

z 

a 

WUCiOOUU 

Ui 

co 

CJ 

M 

u. 

Ui 

UJ 

o 

Vi 

> 

4. 

u  u 

Ui 

X 

Ui  CJ 

z 

Ui 

Ui  u> 

Ui  U) 

Ui 

»- 

Ui 

w  u 

4U 

u. 

Vi  vi 

< 

4 

ViVIViVHJUVI 

Ui 

Ui 

UI 

13 


u 

ac 

< 


CANT  t  IE  VERfcD  dEAP  ANALYSIS  ANO  DESIGN 
INC 


HOoOOOO 


z 

0 tO 

UJmM 
Oh 
ULO 
hU 
zu j 


o»axhO 

4UVKaC3lC 

00 

*“X 

Ui 

<z^>si&a 

h- 

z 

O  N<Z„ 

Z 

h- 

u 

Z  Zxfl-H 

Ui 

3m 

M 

Ui 

z 

CL 

00 

•h 

O 

Ul 

z 

Ul 

Ul  « 

at 

M 

U 

ui  mi 

a.  z 

M 

JWU 

to  Q 

z> 

a 

uiai< 

M 

o 

z 

«ja  * 

Z  K 

Ul 

Mlfito  tUi 

O  < 

at 

ct<  I  aea 

m>N 

Ul 

«/> 

a  »-»-2 

<  Z»u« 

Ui 

O 

-I 

33 

M 

<»- 

< 

40 

► 

Ult- 

W  Ul  3Nr 

(X. 

-l«-> 

-J 

CS  zzz 

*  —  U.I-U- 

Ul 

hO 

M^MHM 

to  MO 

h- 

men 

Z 

•Ulh  (-OL 

•  >Z  Liito 

00 

UM0 

■•JUi-UKQ. 

— 1  JU^^ZUI 

OV4 

UHJQ>Zl 

CJ  <*-*h-eDUJh- 

< 

X 

at ot  *o«z  • 

at  ZhM4i/)< 
h-  X 

h- 

Z 

CJUJUi 

00 

< 

U»>4fMXUUi 

«X 

Z 

< 

uou 

UIO 

UJ 

ZVil  MMtflM 

o 

VJ 

Ui 

oow 

mo 

CO 

mui)u<ua(i*U 

wzui"^«x 

li  (O* 

Oi 

•OO 

•  • 

&OOZ<&<U 

OmjIw>XO 

ZtfHUI  |  HOC 

z'lZt-w.Ofa. 

o 

OQO 

MO 

Q 

m  -*uJ3ax 

UJ 

U13 

a 

CJMPO 

M 

Ul 

uzwviuuui. 

►- 

at  vi 

Ul 

ia  •  * 

at 

aw  uc 

vwtu— a.ujxu.a. 

< 

Ul 

at 

aw 

Ul 

«*WOIOHU4 

X 

X 

Ul 

» 

hOOUCXK 

»- 

M 

UJ 

> 

Ml 

"  II 

Ul 

wa  "a 

« 

K 

MU 

"4 

o«w  at  fit  at® 

md 

00 

— J 

(MW* 

ftJCfittlfiKUhU 

3UI 

Ul 

►- 

M 

hUS«BS33 

U4 

09 

r-jc 

zwtx  taw 

_i«j— »\iPi'rif\'0 

* 

a. 

z 

*”U 

U43333ZIU 

UgZZZZMQ 

u> 

• 

z 

M 

« 

u 

—1 

4auj 

aax 

Nl 

M 

—I 

moe 

oo 

OOO 

z 

* 

< 

OUlUi 

ooo 

oo 

po 

o 

o 

ssss 

o 

V) 

Ui 

ooo 

«o 

o 

0040 

*- 

•J 

OD 

•po 

•  • 

• 

ooo  • 

oou 

(MO 

o 

090010 

V) 

« 

O 

o-*<* 

H 

o 

_i 

u 

UI 

N  •  • 

o 

•  •  t 

ec 

oo 

u> 

ooo 

Ui 

u 

> 

Ml 

«  N 

■ 

NOUN 

o 

lil 

A. 

•J 

mi 

</lQC 

+4 

UM>0 

lib 

o 

«W* 

o 

z 

t— 0&-KM 

U 

«*r 

mJ 

w 

MZUiX. 

ex 

o 

<atu 

auz 

> 

(DMOl 

it 

as 

8 

i 

cu 


a 

ac 

►- 

z 

CO 

CO  COCOCO 

a 

%/i 

CO 

co 

o 

CO 

CO 

CO 

wo 

wo 

V) 

co 

»- 

V) 

VI 

CO 

z 

z 

wo 

VI 

wo 

►-4 

O 

t-4 

to  WO  WO  wo 

co 

ac 

CO 

WO 

a. 

Hi 

UJ 

o 

H- 

UJ 

UJ 

z 

Z3 

LU 

UJ 

< 

X 

o 

n 

co 

Ui 

Ui 

►H 

z 

UJ 

UJ 

Ui 

ac 

h* 

CO 

4 

_i 

UJ 

UJ 

UJ 

Ui 

-J 

UJ 

Ui 

Ui 

a 

Z 

D 

CO 

<1 

UJ1111I UL  immuj 

M 

o 

> 

O 

WO 

K- 

z 

z 

UJ 

»* 

«J  1 

a 

ac 

CO 

Ui 

4 

o 

«J 

■3  | 

H 

a  Q-0.0- 

00 

«a 

< 

4 

b 

a. 

a. 

a 

ac 

-X 

h- 

X 

p-  -»e\iu.  , 

c4 

a 

u. 

H 

4 

<  a 

Ui  l 

M 

a. 

u. 

a 

o 

VI 

a 

z 

c-> 

5 

u. 

u 

** 

X 

1 

a. 

a. 

■J 

u. 

z 

-J 

•J 

4 

t  1 

£ 

a.  0.0. 0.0. 0.0. 

4 

UJ 

Ui 

* 

o 

W 

•-U 

r- 

z 

ac 

aiaj 

UJ 

*4 

»- 

1 

ac 

ac 

z 

OZ 

* 

uoouuoo 

X 

UJ 

UiO 

UJ  ' 

UJ 

o 

Ui 

CJ 

*00«MOKH 

• 

u 

o 

UJ 

UJ^J-J 

-JZ  1 

:•> 

u 

u 

Z 

Ui 

UJ 

Ui 

>ai« 

OJN-  1 

H 

UJ 

c 

UJ 

uj  o 

ac  l 

u 

CJ 

a 

z 

o 

_J4Z 

40. 

UJUJUJUJUJUJUJ 

UJ 

co 

r-»* 

u 

<X 

UJ 

h-«* 

Ui 

UJ 

ZUi 

UI 

Uii 

> 

4 

4 

o 

u 

UJ 

4k 

W«*#» 

•MU 

u 

UJ 

z 

«4 

►4 

(J 

u 

(J 

u 

UJ 

Q 

u 

o 

VJ 

A 

a. 

UL 

1 1 

o 

4 

4 

4 

juuuuuj 

U 

U 

O 

& 

Pv: 


& 


K 


.  * 


5 

►4 

«/) 

Ui 

o 

o 

z 

4 

v> 

•— 

co 

> 


UI 

Ol 

Q 

UI 

(X. 

UI 


NMOOOOO 


n  h  n  h  N  n  n 

c_»><xxi-o 
_jO(/>4  OC-JCD 
<ZZ>  Q.O.  u 
U  <M<ZCl 
Z  ZX&H 


CO  • 

UI  411 
-IcoU 
eOui< 
4-10.  * 

*-*«nco  *ui 

ac<  I  «o 

<h.D<0 

Z^  wi- 
ci 


H>U 


i— ac 

coau>zx 
<XQC  «OmZ  • 

►■Z^htOmO 

U4<_la3c00c*»-0 

*oozza«»g 

4  JUOCLX 
CKZU»tf>UU4Ub»- 
4W  CJZ 

awiuiiU.uuIM 
hUQUOZa 

-14  HO. 

O-'OLQ'CX.tt 

CKOuiujujuif-u; 

l-OSCBBCSS 

Z-IAJC  A*.U.® 

U03D3ZUI 

UUZZXZMU 


UI 

u 

4 

Q. 

to 


u  4 

►h>-M 

o*-*x 


U  CO 

Z  »4  U.I-Q. 

CO  *-<0 

•  UlVf 

-I  -JUI^-JZUI 

o  4— t-cqui|_ 
cl  z»-»-«4«^4 


zuj  ni*-*qc  at  ►-* 

u«JiM>rc 
ZQ»>/>  |  •—« 
Z%*h4UI*0. 
uuiM&uixaa 


JUI 

U3 

44 

<-l> 


a 

a 

uo 


H 

z 

o 

UJ 

a 

Ui 

*— 

—4 

o 

a 

z 

“3 

V 

i 

X 

•— 

03 

o 

CL 

X 

xo 

mU 

OO 

111 

«-»>»■ 

1“  • 

X  • 

CD  • 

• 

> 

4 

UO 

0.0 

<o 

a 

• 

u 

Z 

-1 

N  II 

-J 

*7 

— I 

03 

— 

2 

o 

z 

zo 

o 

H- 

►“ 

»— 

—1 

— 1 

1- 

D 

4 

4 

X 

ls| 

U 

4 

4 

M 

UJ 

A 

y- 

X 

X 

O 

CL 

o 

U'O 

oO 

II 1— 

1-0 

1-  • 

X  • 

_l  • 

>•4 

z 

•— 

uo 

KO 

40 

I— 

Kl 

X 

U 

z 

Ulco 

a 

-1 

■^ui 

u 

4 

z 

03 1- 

o 

o 

u< 

» 

coo 

M 

u 

o 

z 

1- 

U.H 

or 

4 

uu 

UI 

A. 

z 

INI 

z 

z 

at 

«>- 

EX. 

X 

M 

=> 

D 

UI 

UU 

—4 

X 

X 

CL. 

U. 

CBUJ 

M 

a 

oo 

-JO 

EXJO 

Z 

*-.» 

W 

zo 

o  • 

1—  • 

4  • 

*4 

0*4 

o 

uo 

uo 

U<J 

Zr 

to 

— i 

z 

4 

CL 

u 

UJO 

UI 

»-» 

-IUI 

y- 

h- 

OQZ 

Ul 

X 

4 

4w 

X 

4 

IM 

n* 

4 

AW 

►4 

aa 

CL 

U* 

X 

4  UI 

4 

M 

*4 

U. 

K 

4# 

Z 

a 

-itt. 

z 

*- 

O 

u 

4*- 

-4 

Z 

X 

u 

031- 

A 

•— 

uo 

-iO 

uo 

* 

U— 1 

Z 

«— 

o  • 

y-  • 

UI  • 

-O 

a 

a. 

uo 

uo 

oo 

# 

ox 

a 

i— i 

FIG.  3  -  CONT 


Z  trMfWpa 
O  OOOU 


W  LUtULULU 

<QCOOOW 

MQOQOW 

“KOOOO 

JUOOOO 

XU-  •  •  •  • 

0£  UOOU 

o 


ICMAm«N 

UUUO 

UiUMuUJ 

0COOOOO 

1DZOOOO 

azouLio 

0.00000 

_)OJ 

•  •  •  • 

OOOU 


UJ 

Z  sO'O'O'O 

O 

CD 

M 

«»U 

»—  UJLLtUiUJ 

JU 

<00000 

■~o 

NOOOOO 

“hOOOO 

> 

joh^h 

OO  •  •  •  • 

< 

ce  oqoo 

z 

o 

OC<-J 

z 

UJm-4 

4 

Oj-4 

uo 

'U<O<O<0 

KU 

Z^J 

mx 

UJUJU/UJ 

H 

UJ 

acuoooo 

«/> 

UIZOOOO 

\n 

fc 

*30000 

ec 

z 

1- 

UJ 

UJ 

X 

a.  cm 

at  -*«m 

>0CUJIU 

♦  •  •  • 

UJ 

1  oo 

JOOO 

OOOU 

X 

ec 

at-oo 

till 

< 

UJ  UJUJ 

“(JOO 

ac 

uSjZOO 

i-<oo 

-JU-M^ 

< 

u. 

z-ia8§o 

3  •  • 

z  oo 

v>ac 

o 

0«Ud*A(M 

UJUJQOOOO 

UJ 

MM  •  • 

z 

»-x  oo 

M 

M 

*«* 

• 

t—  _J 

XUJ 

-JO 

z 

Z 

oc 

a3 

<z 

a 

M 

¥- 

UJ 

OJ 

M 

«HIN 

wo 

U.4  OH 

O  *«*<N 

J- 

aC4 

z 

Z>  oo 

-JOC 

4 

J-CO  lOOOO 

o 

M 

OX 

X 

W0C30C 

u 

-J  UJUJ 

WO  > 

a 

z-*< 

ui<aaoo 

UJ 

o 

OO-' 

u. 

UMUIZOO 

-J 

u. 

u 

a 

Oh*3UO 

CO 

z 

<mUOUO 

4 

M 

4-J-« 

ac 

•«z-»«h»vx 

M  • 

«c 

UJ 

at—  •  • 

ee>  • 

J* 

CU 

UJ 

4  oo 

z 

UJO  CL 

X 

>o 

>  >z 

M 

«— J< 

z> 

oc 

a 

4 

40 -> 

z 

zuj 

z 

a 

UiK*  « 

UJ 

►- 

UJ 

— J 

*- 1  > 

M 

wo 

ec 

4 

VtC  MN 

V) 

z 

UJ  a  ■■  CM  CM  4 

UJQ  «0 

UJ  OHM 

w 

x 

u 

uiuz 

U  M 

a 

UJ 

z 

IX 

M 

M 

o 

o 

UJ 

oc 

UJ 

UJ 

-J 

O 

a) 

< 

«u 

ac 

uu 

u 

hU 

H- 

<«■> 

w 

> 

WTNwu 

< 

OU 

4 

J- 

X 

OUJUJ 

oo 

4 

z 

< 

ooo 

oo 

u 

UJ 

OOO 

wr\o 

_Jm(7> 

OJ 

•oo 

•  • 

WJ 

<W>— "* 

ooo 

CMO 

UJ 

UO-* 

UJ 

0-*«”l 

H 

M 

acu 

IM 

in  •  t 

■*4. 

UJ 

oc 

oo 

X 

X 

UJ 

M 

UJ 

> 

N  H  M 

H  M 

f- 

IU 

VJ 

-j 

UJ 

*- 

■3f» 

H* 

* 

am 

z 

-J 

« 

M 

UJ 

xcikUi 

UJX 

63 


FIG.  3  -  CONT 


CAKTILEVEREC  BEAM 


moo 

infn*-* 

oo 

UQU 

Oujuj 

oo 

u 

UJUIUJO 

OUO 

oo 

w 

o©r**o 

OUW 

inu 

w 

uo'uo 

•uu 

•  • 

• 

UO«  • 

OUO 

Olu 

u 

auocg'T 

0-4JW 

(J 

«r^>sr 

"1  •  I 

u 

•  •  • 

wo 

m 

uvo 

R  R 

w 

™  H  ™  H 

toot 

xi 

mi 

<o.uj 

ecx 

O 

euTlaz 

*  * 

*  * 


*  * 


* 

Z 

* 

w 

• 

►4 

* 

4- 

* 

< 

* 

Kl 

« 

M 

* 

QL 

X 

* 

o 

fr* 

* 

<L 

z 

* 

* 

A 

X 

X 

« 

« 

* 

M 

a 

z 

« 

o 

n 

* 

X 

o 

4-4 

* 

<r 

»- 

* 

X 

u. 

v-A 

* 

z 

* 

og 

o 

z 

o 

o 

* 

u. 

* 

u 

at 

— UJ 

* 

t- 

a 

« 

7 

aro 

a 

UJ 

O 

OO 

* 

u 

z 

■* 

M 

4-0 

u. 

*-A 

»- 

oo 

* 

«* 

♦ 

u- 

X 

• 

♦ 

►- 

IT 

at 

1  o 

n 

X 

<r 

♦ 

LL 

mm 

V-J 

z 

o 

o 

♦ 

M 

«/IU 

UJ 

UJ 

* 

# 

Z 

a 

-JUJ 

O 

o 

mo 

♦ 

W 

** 

o 

*»o 

B- 

o 

— o 

* 

♦ 

o 

o 

0.10 

X 

\r\ 

<1  fM 

♦ 

Z3 

• 

• 

a 

o 

o 

*  4 

4*4 

1*4 

4*4 

i  * 

z 

—I 

W 

-■a 

M 

H 

VI 

4- 

“} 

o-* 

Z 

CL 

UJ 

c 

(_> 

CCIS1RAINT  VALUES  IC-VECTCR) 

II  C  « 14CCCE  01  -C.940C0E  00  0.32f*7E  01  -0.600CCE  00 


com 

22 

OO 

» ii»  1 

>«■ 

l-(- 

0 

• 

a cot 

Ul 

ujuj 

CO 

Y-y- 

— *« 

SJ> 

<0 

t 

0 

aroc 

1 

00 

U.U. 

Ef  . 


S.T; 

•1 


s 


I 


4 


O 

Z2 

1 

33 

Ui 

U.“* 

>r 

_iCO 

in 

UiX 

IM 

DU 

rr> 

in 

ni 

K 

zz 

• 

Z 

0 

l 

to 

M 

< 

VI 

z 

OC 

<V 

0 

►— 

v— 

K 

!/)</> 

a 

u 

z 

■«* 

cn 

COCCI 

M 

a c 

Z 

Ui1" 

z 

•*"11 1 

Ui 

Y- 

O 

_l-i 

0 

—-u 

oc 

in 

O 

M 

Of"* 

0C<O 

Y- 

Z 

1- 

o«r 

«/> 

O 

Ui 

0m 

«* 

Ooo 

►-«n 

z 

u 

O 

■  '■»- 

X 

wm 

US' 

u 

•— 

Vi  —1 

u. 

?  • 

ui  • 

Ui 

u 

»i 

-J  — 

a 

«  O 

>0 

Ui 

Ui 

CD  1 

u. 

X 

1  1 

war 

1- 

Ui 

zn- 

Z 

0 

< 

> 

CJV- 

*■— 

c 

— < 

-J 

— 

—*-» 

tn«J 

O 

»-«n 

0 

l~ 

ac-ucr> 

Z 

Ui 

UI 

in  1 

uoc 

M 

0 

Ui  I  0 

O 

ao 

OUI 

ujuj 

«Ui 

> 

1 

— V 

<** 

ujw> 

m 

V- 

u* 

«ti* 

—>'0 

IVI 

0 

0 

QC^-U 

O 

l-t-i 

3 

OCOu> 

IM 

-o 

IX  CO 

Z 

O'* 

— 

<0 

«*■«* 

zT  ai 

X 

• 

>  • 

*-  • 

C 

u-u 

»"l 

w 

0 

Z<3 

wz 

ui 

Ui 

—  II—  VM 

*- 

z 

"*  1 

at  i-4 

oc 

OC 

»-inv) 

a 

0 

< 

« 

n 

«ai<ii 

O 

Mt 

ae 

at 

z««c 

H 

«/> 

t- 

UiH« 

Ui 

IV 

— 

V)« 

atm 

ae 

OC 

X 

“9 

0** 

Z-* 

UiZ 

Ui 

Ui 

oc 

z 

aj 

UI 

O 

XU 

X 

X 

Ui 

** 

O 

O 

O 

1-0 

V- 

Y- 

h- 

65 


NLMeiR  OF  ITERATICF5  *  9 

OBJECTIVE  FUNCTION  NAS  EVALUATED  35  TIKES 

CONSTRAINT  FUNCTIONS  NERE  EVALUATED  35  TIMES 


I 

» 


7  w*v*  1.1  vw;  r.  a  a:  -T^*r 


t: 

L’ 


i 

A 

A 

A 

,N 

A 

A 


B 


oo 


UPPER 

BOUNC 

0000E 

OCCOE 

U1«\l 

•  * 

oo 

OOOO 

a 

QCZULIiiiUJuj 

uj-iuoua 

d|IN 

0.00000 

OO 

asuooo 
3  OOOO 

04 

UlrnUJ 

_J0>0> 

•  •  •  • 

J|*>N 

OOOO 

UJ 

■* 

>CD<X> 

O' 

a 

•  • 

IOCOOpH 

•0 

OO 

OOOO 

UI 

•  • 

OUJUJUJUJ 

o 

o~« 

jmoiPN 

zr-on^-o' 

UJ 

oo 

>crn'u  <r 

o' in  r*  O' 

-1 

q-QUIUJ 

M^-O'O' 

< 

> 

UJZOO 

*300 

OOOO 

oood 

Z 

JCOOU 

o 

Ifl"* 

•  • 

««<o<o 

H 

oo 

•-4.-IPU-U 

u 

a 

z 

arZujujUJuj 

3 

UJ30000 

u. 

t 

JCCJOOOO 

-JO 

0(00000 

«z 

1  P4HM«H 

(U  -J'X 

*Mt>«fa4PU 

u  • 

•  •  •  • 

—lac 

OOOO 

VI 

cn 

U<1 

•  III 

►- 

-* 

z 

V) 

• 

:d 

o 

H* 

a 

s/l 

M 

«/) 

-JZ 

UJ 

KZ 

UJ 

< 

at 

uo 

— j 

• 

n 

co  »'*t/>'Or* 

Zl-* 

CO 

>  • 

a 

aac 

z 

<« 

0-«0l 

-jz 

u 

U.Z 

•z 

VI 

++ 

o 

ar 

u 

H- 

oju 

■«* 

<-/ 

«. 

;>-J 

> 

u 

M 

M 

1--* 

z 

<C 

u 

at 

o< 

o 

Q-<fM 

W 

H-Mmf 

** 

UJCO 

*u 

•M 

M 

r-* 

*X7 

w 

V# 

a 

OJ 

UJ 

Ui 

u 

DO 

a 

fcj 

in  oo 

uvoo 

oo 

ooo 

X 

OUJUJ 

COD' 

in 

UJUJUJO 

z 

QOQ 

r~ 

UVJlflU 

UJ 

ooo 

tVf«4 

cO 

M/'P'O 

tt) 

•oo 

•  • 

• 

• 

OoO 

HOD 

0* 

tj'inr-o 

a 

o-jfn 

H 

o 

— <*r  cr  — < 

UJ 

*  • 

<o 

*  •  » 

a 

oo 

•o 

OOO 

UJ 

> 

DIM 

N  If 

M 

It  It  H  N 

UJ 

UI* 

UJ*-< 

K 

OLW*t- 

Z 

— J 

c-a-JCD 

Z 

u 

V»lujv 

u 

<CLUJ 

cux 

> 

CDV'QX 

S5 

o 

o 

I 

m 


O 


66 


PRCGPAP  CALLS  TO  ANAL  IZ 
IC/LC  CALLS 


67 
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